BROAD-BEAM IMAGING METHODS 



CROSS-REFERENCES TO RELATED APPLICATIONS 
[0001] This application is a continuation— and claims the priority benefit— of 
5 commonly owned U.S. patent application number 10/211,391 entitled "Broad-Beam 
Imaging/' filed August 1, 2002 which, in turn, claims the priority benefit of commonly 
owned U.S. provisional patent application number 60/370,632 entitled "Broad-Beam 
Imaging," filed April 5, 2002. The subject matter of these two applications is 
incorporated herein by reference. 

10 [0002] Commonly owned U.S. patent application number 10/211,391 entitled "Broad- 
Beam Imaging" is also a continuation-in-part— and claims the priority benefit— of 
commonly owned U.S. patent application number 10/039,922 entitled "Block Switching 
in Ultrasound Imaging," filed October 20, 2001. The subject matter of this application is 
incorporated herein by reference. 

15 [0003] Commonly owned U.S. patent application number 10/211,391 entitled "Broad- 
Beam Imaging" is also a continuation-in-part— and claims the priority benefit— of 
commonly owned U.S. patent application number 10/039,862 entitled "Simultaneous 
Multi-Mode and Multi-Band Ultrasonic Imaging," filed October 20, 2001. The subject 
matter of this application is incorporated herein by reference. 

20 [0004] Commonly owned U.S. patent application number 10/211,391 entitled "Broad- 
Beam Imaging" is also a continuation-in-part— and claims the priority benefit— of 
commonly owned U.S. patent application number 10/039,910 entitled "System and 



Method for Coupling Ultrasound Generating Elements to Circuitry," filed October 20, 
2001 which, in turn, is a continuation-in-part of U.S. patent application number 
09/860,209, now U.S. patent number 6,569,102, entitled "Miniaturized Ultrasound 
Apparatus and Method," filed on May 18, 2001 which, in turn, is a continuation of U.S. 
5 patent application number 09/378,175, now U.S. patent number 6,251,073, entitled 
"Miniaturized Ultrasound Apparatus and Method," filed on August 20, 1999. The 
subject matter of these applications is incorporated herein by reference. 
[0005] This application is related to co-pending U.S. patent application number 
09/872,541 entitled "System and Method for Phase Inversion Ultrasonic Imaging," filed 
10 May 31, 2001. The subject matter of this commonly owned application is hereby 
incorporated by reference. 

[0006] This application is further related to co-pending U.S. patent application 
number 10/101,661, now U.S. patent number 6,663,567, entitled "System and Method for 
Post-Processing Ultrasound Color Doppler Imaging," filed March 19, 2002. The subject 
15 matter of this commonly owned application is hereby incorporated by reference. 

BACKGROUND 

FIELD OF THE INVENTION 

[0007] The invention is in the field of imaging and more specifically in the field of 
20 ultrasonic imaging. 
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PRIOR ART 

[0008] Ultrasonic imaging is a method of analysis used for examining a wide range of 
materials. The method is especially common in medicine because of its relatively non- 
invasive nature, low cost, and fast response times. Typically, ultrasonic imaging is 
5 accomplished by generating and directing an ultrasound beam into a material under 
investigation in a transmit phase and observing reflections generated at the boundaries 
of dissimilar materials in a receive phase. For example, in medical applications 
observed reflections are generated at boundaries between a patient's tissues. The 
observed reflections are converted to electrical signals (channel data) by receiving 
10 devices (transducers) and processed, using methods known in the art, to determine the 
locations of echo sources. The resulting data is displayed using a display device such as 
a monitor. 

[0009] The prior art processes of producing an ultrasound beam and analyzing 
resulting echoes is called "beam forming." The production process optionally includes 

15 defining "transmit" beam characteristics through aperture apodization, steering, and/or 
focusing. The analysis process optionally includes calculating a "receive beam" wherein 
received echoes are processed to isolate those echoes generated along a narrow region. 
This calculation includes the identifying one-dimensional line along which echoes are 
assumed to have been generated, and is therefore referred to herein as "echo line 

20 calculation." Through beam forming a one-dimensional set of echolocation data is 

generated using each transmit and/or receive beam. Echolocation data is positional data 
relating to the physical location of one or more echo source and optionally includes 
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intensity, velocity and/or similar physical information. Echolocation data may include 
post-beam forming raw data, detected data, or image data. Multidimensional 
echolocation data, such as an ultrasound image, is generated by scanning a field of view 
within the material under investigation using multiple transmit and/or receive beams. 
5 [0010] The ultrasound beam transmitted into the material under investigation during 
the transmit phase is generated by applying electronic signals to a transducer. The 
ultrasound beam may be scattered, resonated, attenuated, and/or reflected as it 
propagates through the material under investigation. A portion of the reflected signals 
are received at transducers and detected as echoes. The receiving transducers convert 
10 the echo signals to electronic signals and optionally furnish them to an echo line 

calculator (beam former) that performs the echo line calculation inherent to analysis 
using a receive beam. 

[0011] After beam forming, an image scan converter uses the calculated echolocation 
data to generate image data. In prior art systems the image formation rate (the frame 

15 rate) is limited by at least the total pulse return times of all ultrasound beams used to 
generate each image. The pulse return time is the time between the transmission of the 
ultrasound beam into the material under investigation and the detection of the last 
resulting reflected echoes. The limited frame rate may result in temporal artifacts 
caused by relative movement between the ultrasound system and a . material under 

20 investigation. 

[0012] FIG. 1 shows a prior art ultrasound system, generally designated 100. 
Ultrasound system 100 includes an element array 105 of transducer elements 110, a 



backing material 120, an optional matching layer 130, a transmit/receive switch 140 and 
a beam transmitter 150. Backing material 120 is designed to support element array 105 
and dampen any ultrasound energy that propagates toward backing material 120. 
Matching layer 130 transfers ultrasound energy from transducer elements 110 into the 
5 material under investigation (not shown). Transducer elements 110, include individual 
transducer elements 110A-110H individually coupled by conductors 115 and 117, 
through transmit/receive switch 140, to a beam transmitter 150. Transmit/receive switch 
140 may include a multiplexer 145 that allows the number of conductors 117 to be 
smaller than the number of conductors 115. In the transmit phase, beam transmitter 150 
10 generates electronic pulses that are coupled through transmit/receive switch 140, 

applied to some or all of transducer elements 110A-110H, and converted to ultrasound 
pulses 160. Taken together, ultrasound pulses 160 form an ultrasound beam 170 that 
probes the material under investigation. 

[0013] Ultrasound beam 170 may be focused to limit the region in which echoes are 
15 generated. When echo sources are restricted to a narrow region the calculation of echo 
location data may be simplified by assuming that the echo sources lie along a " transmit 
line/ 7 With this assumption, the task of the echo beam calculator is reduced to a 
problem of determining the position of an echo source in one dimension. This position 
is established using the return time of the echo. The accuracy of this assumption and the 
20 spacing of transmit lines are significant factors in determining the resolution of prior art 
ultrasound systems. Finely focused beams facilitate higher resolution than poorly 
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focused beams. Analogous assumptions and consequences are found in analyses 
involving calculated receive beams. 

[0014] FIG. 2 shows a prior art focusing system in which element array 105 is a phased 
array configured to focus ultrasound beam 170 by varying the timing of electronic 
5 pulses 210 applied to transducer elements 110A-110H. In this system, electronic pulses 
210, are generated at beam transmitter 150 and passed through transmit/receive switch 
140. Electronic pulses 210 are delayed using a delay generator (not shown) and coupled 
to transducer elements 110A-H. Ultrasound beam 170 is formed when transducer 
elements 110A-H convert properly delayed electronic pulses 210 to ultrasound pulses 

10 160 (FIG. 1). Once formed, ultrasound beam 170 is directed along a transmit beam line 
250 including a focal point 230 with a resulting beam waist 240 characterized by a width 
of ultrasound beam 170. In a similar manner phased excitation of element array 105 is 
used to direct (steer) ultrasound beam 170 in specific directions. The cross-sectional 
intensity of ultrasound beam 170 is typically Gaussian around a focal point and includes 

15 a maximum along transmit beam line 250. The shape of ultrasound beam 170 may 
depend on aperture apodization. 

[0015] In a scanning process, ultrasound system 100 sends a series of distinct 
ultrasound beam 170 along another, different transmit beam line 250 to form an image 
over more than one spatial dimension. A specific ultrasound beam 170 is optionally 
20 transmitted in several transmit/receive cycles before generating another ultrasound 

beam 170. Between each transmit phase a receive phase occurs, during which echoes are 
detected. Since each ultrasound beam 170, included in an ultrasound scan, requires at 
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least one transmit/receive cycle the scanning processes may take many times the pulse 
return time. This pulse return time, determined by the speed of sound in the material 
under investigation, is a primary limitation on the rate at which prior art ultrasound 
images can be generated. In addition, undesirable temporal anomalies can be generated 
5 if transducer elements 110A-110H move relative to the material under investigation 
during the scanning process. 

[0016] FIGs. 3A through 3E show a prior art scanning process in a phased array 310 of 
eight transducer elements, designated 110A through 110H. Subsets 320A-320E of the 
eight transducer elements 100A-110H are each used to generate one of distinct 

10 ultrasound beams 170A-170E. For example, FIG. 3A shows ultrasound beam 170A 
formed by subset 320A, including transducer elements 110A - HOD. The next step in 
the scanning process includes forming ultrasound beam 170B using subset 320B 
including transducer elements HOB - HOE as shown in FIG. 3B. In this example, a 
transmit beam line 250B associated with ultrasound beam 170B passes through a focal 

15 point 230B, which is displaced from a focal point 230 A by a distance typically equal to 
the width of one transducer element 110. As shown by FIGs. 3C through 3E, each subset 
320C through 320E, used to produce each ultrasound beam 170C through 170E, is 
displaced by one transducer element 110 relative to subsets 320B through 320D, 
respectively. Echoes detected in the receive phase, occurring between each transmit 

20 phase, are used to generate echolocation data and these echolocation data are typically 
combined to form an image suitable for display. The scan process may be repeated to 
produce multiple images. 



[0017] In practice, phased array 310 may include sixty-four, one hundred and twenty- 
eight, or more transducer elements 110. The resolution of the echolocation data depends 
on the aperture and the number of transducer element 110, and on the degree to which 
transmit beam line 250 accurately represents possible echo sources within ultrasound 
5 beam 170. Representation of ultrasound beam 170A-E using beam line 250A-E is an 
approximation that determines the resolution of resulting echolocation data. A poor 
approximation will limit the resolution of the resulting echolocation data. A maximum 
width of ultrasound beam 170A-E is, therefore, limited by the desired resolution of the 
echolocation data. The accuracy of the approximation is a function of distance from 

10 focal points 230A-E, the approximation being less accurate at further distances. 

[0018] Common practice includes generating several ultrasound beams with different 
focal point 230 A-E, and using each set of received echoes to generate data near focal 
points 230A-E. Prior art data generation may be limited to an area near focal points 
230A-E because, at further distances, the transmit beam line 250 approximation may not 

15 be sufficiently accurate to provide the echolocation data of a desired resolution. 

Typically one receive or transmit beam line 250 is generated for each transmit/receive 
cycle. The number of beams required to image an area is dependent on both the width 
and depth of the area to be imaged as well as the desired resolution. By using only 
echoes near focal point 230, only a small portion (e.g. < 10%) of the total received signal 

20 is used, with the remainder of the received signal being discarded. The prior art makes 
inefficient use of detected signal. Similar disadvantages occur in systems utilizing 
synthetic receive lines. 



[0019] In the prior art the area to be covered, transmit beam width, number of 
transmit beam 170, and echolocation data resolution are interdependent. The transmit 
beam width determines the minimum lateral resolution width of the echolocation data. 
Since each transmit beam 170 covers only a limited area, a greater number of transmit 
beam 170 are required to image a larger area. Use of a greater number of transmit beam 
170 lengthens the minimum time required to generate an image. 
[0020] Disadvantages of the prior art, such as an image formation rate restricted by 
pulse return time and inefficient signal use, have prevented prior art ultrasound systems 
from taking full advantage of advances in micro-processing power. The prior art 
endures these disadvantages in order to generate images with the highest possible 
resolution. 
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BRIEF DESCRIPTION OF THE DRAWING 
[0021] FIG. 1 shows a prior art ultrasound system; 
[0022] FIG. 2 shows a prior art method of focusing an ultrasound beam; 
[0023] FIGs. 3A through 3E show a prior art scanning process using a phased array of 
5 eight transducer elements; 

[0024] FIG. 4 is a flow chart showing an overview of a broad-beam method according 
to an embodiment of the invention; 

[0025] FIG. 5 shows a block diagram of a broad-beam system according to an 
embodiment of the invention; 
10 [0026] FIG. 6 is a flow chart showing details of a broad-beam design step according to 
an embodiment of the invention; 

[0027] FIG. 7 A shows an ultrasound beam generated using a linear transducer array 
according to an embodiment of the invention; 

[0028] FIG. 7B shows an ultrasound beam generated using a curvilinear transducer 
15 array according to an embodiment of the invention; 

[0029] FIG. 7C shows an ultrasound beam that results in an insonified region 
generated according to an embodiment of the invention; 

[0030] FIG. 7D shows a plot of ultrasound intensity through a cross-section of a 
broad-beam; 

20 [0031] FIG. 8 is a flowchart showing details of a transmit step according to an 
embodiment of the invention; 
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[0032] FIG. 9 is a flowchart showing details of a receive step according to an 
embodiment of the invention; 

[0033] FIG. 10 shows stored data arranged in a channel data array according to an 
embodiment of the invention; 
5 [0034] FIG. 11A shows an echolocation data array including a first axis indicating X 
position and a second axis indicating Y position according to an embodiment of the 
invention; 

[0035] FIG. 11B shows an alternative embodiment of the echolocation data array 
including first axis indicating angle (9) and second axis indicating radius (R) according 
10 to an embodiment of the invention; 

[0036] FIG. 12A shows a Cartesian coordinates system including, for the purposes of 
illustration, eleven "X" divisions separating data bins according to an embodiment of 
the invention; 

[0037] FIG. 12B shows a radial coordinate system representing the area insonified by 
15 an ultrasound beam according to an embodiment of the invention; 

[0038] FIGs. 13A and 13B show ultrasound propagating from transducer elements to 
objects within a material under investigation according to an embodiment of the 
* invention; 

[0039] FIG. 14 shows channel data produced from echoes according to an 
20 embodiment of the invention; 

[0040] FIG. 15 shows echolocation data generated using the channel data shown in 
FIG. 14 according to an embodiment of the invention; 
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[0041] FIG. 16 is a flowchart showing a method included in an echo area calculation 
according to an embodiment of the invention; 

[0042] FIG. 17 shows a graph illustrating three alternative apodization functions 
according to an embodiment of the invention; 
5 [0043] FIG. 18 shows ultrasound transmitted from Jwo transducer elements and 
striking an ultrasound reflective object; 

[0044] FIG. 19 shows signals generated by an SCE transducer element stored in a 
channel data array according to an embodiment of the invention; and 
[0045] FIG. 20 is a flowchart showing details of an echo area calculation step 
10 according to an embodiment of the invention. 
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SUMMARY OF THE INVENTION 
[0046] Some embodiments of the invention include a method of probing a material 
under investigation comprising the steps of (1) using a plurality of transducers to 
transmit an ultrasound beam into the material under investigation, the ultrasound beam 
5 including components generated by each transducer in the plurality of transducers, (2) 
receiving echoes generated by interactions between the ultrasound beam and the 
material under investigation, (3) generating first data from the received echoes, the first 
data having values that include phase and magnitude information and being 
associatable with a time dimension and distributed over at least one spatial dimension, 

10 (4) using the phase and/or magnitude information to distinguish echoes, among the 

received echoes, resulting from ultrasound beam components generated by a subset of 
transducers in the plurality of transducers, and (5) transforming the first data into 
second data using the distinguished echoes, the second data having values distributed 
over at least one more spatial dimension than the first data. 

15 [0047] Some embodiments of the invention include a method of probing a material 
under investigation comprising the steps of (1) transmitting one ultrasound beam into 
the material under investigation, (2) receiving echoes generated by interactions between 
the ultrasound beam and the material under investigation, (3) generating first data from 
the received echoes, the first data having a value that includes phase and magnitude 

20 information and associatable with time and at least a first spatial dimension, and (4) 
transforming a portion of the first data into second data using a transform capable of 
producing second data distributed over at least a second spatial dimension and a third 
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spatial dimension, the transform using the phase and/or magnitude information to select 
the portion of first data to be transformed. 

[0048] Some embodiments of the invention include a method of probing a material 
under investigation comprising the steps of (1) transmitting one or more ultrasound 
5 beam into the material under investigation, (2) receiving first echoes generated by 
interactions between one of the transmitted one or more ultrasound beam and the 
material under investigation, the interactions occurring at points distributed over at least 
a first spatial dimension and a second spatial dimension, (3) generating first data from 
the received first echoes, the first data having values distributed in a time dimension 

10 and additionally distributed over at least the first or the second spatial dimension, (4) 
transforming the first data into second data having values distributed over at least both 
the first and the second spatial dimension, (5) transmitting an other ultrasound beam 
into the material under investigation, (6) receiving further echoes generated using the 
other ultrasound beam, (7) generating third data using the received further echoes, the 

15 third data being echolocation data and having a dimensionality, and (8) combining the 
third data with the second data, the combination having the same dimensionality as the 
third data. 

[0049] Some embodiments of the invention include a method of probing a material 
under investigation comprising the steps of (1) transmitting an ultrasound beam into the 
20 material under investigation, (2) receiving echoes generated by interactions between the 
transmitted ultrasound beam and the material under investigation, (3) generating first 
data using the received echoes, the first data having values associatable with time and a 
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number of positions in a first spatial dimension, the number of positions being at least 64 
and the association with the number of positions being independent of the association 
with time, and (4) transforming the first data into second data having values associatable 
with at least the first spatial dimension and a second spatial dimension. 
5 [0050] Some embodiments of the invention include a method of probing a material 
under investigation comprising the steps of (1) using a plurality of ultrasound 
transducers to transmit an ultrasound beam into the material under investigation, (2) 
receiving echoes generated by interactions between the ultrasound beam and the 
material under investigation, (3) generating first data from the received echoes, the first 

10 data having a plurality of values associatable with time and with a first number of 

positions in a first spatial dimension, the first number of positions being more than one, 
and (4) generating second data from the first data, the second data having values 
associatable with a second spatial dimension and with a second number of positions in 
the first spatial dimension, the first number of positions being fewer than the second 

15 number of positions; wherein at least one of the values of the second data, associatable 
with one of the second number of positions but not with any of the first number of 
positions, is generated without interpolation between values of the first data. 
[0051] Some embodiments of the invention include method of probing a material 
under investigation comprising the steps of (1) using a plurality of transducers to 

20 transmit a plurality of ultrasound beams into the material under investigation, (2) 

receiving first echoes generated by interactions between a first member of the plurality 
of ultrasound beams and the material under investigation, (3) generating first echo data 
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from the received first echoes, the first echo data having values associatable with a 
temporal dimension and separately associatable with more than one position in at least a 
first spatial dimension, the values including phase and magnitude information, (4) 
receiving second echoes generated by interactions between at least a second member of 
the plurality of ultrasound beams and the material under investigation, (5) generating 
second echo data from the received second echoes, the second echo data having values 
associatable with a temporal dimension and separately associatable with more than one 
position in at least a second spatial dimension, (6) generating first echolocation data 
using the first echo data and a data transform responsive to the phase and/or magnitude 
information, (7) using the second echo data to generate second echolocation data, and (8) 
combining the first and the second echolocation data to produce third echolocation data 
having the same dimensionality as the first echolocation data. 
[0052] Some embodiments of the invention include a method of generating 
echolocation data comprising the steps of (1) generating first data by converting echoes 
into electronic signals, the first data having a plurality of values associatable with time 
and separately associatable with a plurality of positions in at least one spatial dimension, 
the plurality of values including phase and magnitude information, and (2) generating 
the echolocation data using the first data and a data transform responsive to the phase 
and/or magnitude information, the echolocation data having at least one value derived 
from two or more members of the plurality of values associatable with different 
positions in the plurality of positions. 
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DISCLOSURE OF THE INVENTION 
[0053] New broad-beam technologies are systems and methods that allow 
multidimensional (area or volume) echolocation data to be generated from as few as one 
5 ultrasound beam. These technologies include generating an ultrasound beam and 
transmitting it into a material under investigation, generating echo signals from 
resulting echoes, and processing the echo signals to produce echolocation data 
distributed in two or more dimensions. 

[0054] Broad-beam technologies are less complex than prior art ultrasound systems 
10 and methods. For example, broad-beam systems and methods are not restricted by the 
use of transmit lines, scan lines or receive lines, and broad-beam systems and methods 
can generate multidimensional echolocation data from as few as one transmitted 
ultrasound beam. Dependence on transmit lines and receive lines is eliminated because 
broad-beam technologies do not require an assumption that echo sources are located 
15 along a one-dimensional line, such as transmit beam line 250 and/or a receive line. 

Broad-beam systems and methods do not require multiple beam scanning or scan lines 
to generate a two dimensional image. Also, unlike the prior art, the resulting 
echolocation data may result from a single transmitted ultrasound beam that may be 
distributed over two dimensions. Using broad-beam systems and methods, a majority 
20 of the received echo signals may be used for image generation. 

[0055] Unlike prior art embodiments, broad-beam systems and methods do not 
necessarily depend on a transmitted ultrasound beam's shape or width to determine the 
resolution of echolocation data. This independence arises because broad-beam systems 
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include no assumption that a transmitted ultrasound beam is approximated by a 
transmit line or a column surrounding a transmit line. Generally, ultrasound beams 
(broad-beams) used in broad-beam systems and methods are wider than the finely 
focused ultrasound beam 170 used in the prior art. 
5 [0056] Broad-beam systems and methods manipulate data differently than the prior 
art. Broad-beam systems and methods are based on multidimensional de-convolution 
algorithms that convert echoes received at receiving transducers into echolocation data, 
thereby generating multidimensional echolocation data from a single transmitted 
ultrasound beam. For example, in one embodiment a de-convolution algorithm 

10 (calculation) affects a transform from two dimensional (time, ultrasound transducer) 
raw data to two dimensional (X,Y position) echolocation data. The two dimensional 
(time, ultrasound transducer) raw data is optionally generated using a single 
transmitted ultrasound beam, and without assuming a transmit line or a receive line. 
The two dimensional echolocation data is distributed over an area requiring at least two 

15 spatial dimensions for representation. The data manipulation included in broad-beam 
systems and methods is capable of using a single transmitted ultrasound beam to 
produce a two-dimensional image configured for display on a display device. 
[0057] Broad-beam systems and methods take advantage of increases in micro- 
processor power and advances in integrated circuit technologies. Current micro- 

20 processors are capable of performing broad-beam data analysis at a rate that is faster 
than the rate at which individual ultrasound beams can be transmitted and received 
using prior art beam-forming technologies. While prior art technologies are restricted 
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by the pulse return time and the number of individual ultrasound beams needed to 
image an area, embodiments of the broad-beam approach leverage ongoing advances in 
computing technology. Broad-beam systems and methods achieve image generation 
rates that are not primarily limited by the use of narrowly focused ultrasound beams, as 
5 in the prior art. 

[0058] For example, in a conventional system imaging to a depth of 200 mm, 128 
transmit/receive cycles require 33.3 milliseconds based on a speed of sound of 1.54 
mm/microsecond. This rate yields a frame rate of approximately 30 frames/second with 
an image resolution across the image area, perpendicular to the axis of element array 

10 105, of 128 lines. In comparison, using an embodiment of the invention to image the 
same depth, a similar resolution can be obtained using five to seven transmit/receive 
cycles requiring a total of 1.3 to 1.8 milliseconds. These times limit the resulting frame 
rate to 769 and 549 frames/second respectively. In various embodiments, images, with 
image resolutions of 128 lines as above, are obtained in less than 25, 17, 10, 5, or 2 

15 milliseconds. 

[0059] Some embodiments of broad-beam technologies result in images that minimize 
the occurrence of undesirable temporal anomalies associated with prior art scanning 
processes. The multidimensional echolocation data derived from a broad-beam 
ultrasound beam is representative of a section of the material under investigation during 
20 the short period of a pulse return time. Since this time is shorter then the time required 
to accomplish a two-dimensional (multiple beam) scan in the prior art, the probability of 
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relative movement between the transducers and the material under investigation during 
the data collection is reduced relative to the prior art. 
[0060] Broad-beam systems and methods do not depend on the prior art 
approximation that an ultrasound beam can be represented by a line, such as beam line 
5 250. Therefore, the resolution of resulting echolocation data is not a function of distance 
from a focal point, such as prior art focal points 230A-E. Broad-beams are typically 
wider, and capable of imaging areas larger, than each of the focused beams of the prior 
art. 

[0061] Since each broad-beam is capable of imaging an area larger than prior art 
10 ultrasound beams, the number of ultrasound beams required to image a specific area is 
reduced relative to prior art. Because fewer, such as only one, ultrasound beams are 
required, broad-beam systems and methods may use less power to image a material 
under investigation than prior systems. Using less power decreases the amount of 
energy deposited in the material under investigation, and decreases the amount of 
15 electricity required to generate each image. Reduced electrical requirements may benefit 
devices using self-contained power sources, such as batteries. 

[0062] Embodiments of broad-beam technology include an area forming™ process of 
producing, receiving, and analyzing an ultrasound beam wherein a set of echolocation 
data, distributed over an area requiring two spatial dimensions for representation, is 
20 generated using as few as one ultrasound beam. The receive points at which echo 

detection occurs and echolocation data is generated may be anywhere within the probed 
region. The receive points optionally lie along a variable grid whose granularity and 
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regularity vary with position. Other embodiments of broad-beam technology include a 
volume forming process similar to area forming except that three spatial dimensions are 
required to adequately represent the echolocation data generated using as few as one 
ultrasound beam. Area forming and volume forming are optionally combined with 
5 non-spatial dimensions, such as time and velocity to achieve multidimensional forming 
processes. 

[0063] FIG. 4 is a flow chart showing an overview of a broad-beam method according 
to an embodiment of the invention and generally designated 400. Method 400 begins 
with a broad-beam design step 410 that includes determination of the number and 

10 shapes of ultrasound beams (broad-beams) needed to image an area or volume. Within 
this step, desired characteristics of at least one of the determined broad-beams are 
calculated and parameters for the broad-beam's generation are established. The desired 
characteristics of each broad-beam may include factors such as position, direction, 
width, intensity, dispersion, or the like. The parameters may include voltages, aperture 

15 functions, excitation delays, and such. 

[0064] In a transmit step 420, the broad-beam designed in step 410 is generated and 
transmitted into a material under investigation. Transmit step 420 includes generation 
of an electronic waveform using, for example, a digital or analog waveform generator. 
This waveform is coupled to multiple channels, each of which may be independently 

20 delayed and amplified using devices such as a multi-channel delay generator and a 
multi-channel power amplifier. Typically, delay times are selected responsive to the 
desired shape, width and direction of the broad-beam. The amplified waveforms excite 



transducer elements 110 causing the broad-beam to be transmitted into a material under 
investigation. 

[0065] A receive step 430 uses transducer elements 110 to detect echoes produced by 
the transmitted broad-beam. Transducer elements 110 generate electronic signals 
5 responsive to the detected echoes. The resulting electronic signals (analog channel data) 
are optionally filtered using an analog filter and digitized, typically with a multi-channel 
A/D converter, to generate digital channel data. In one embodiment, the channel data 
preferably includes both amplitude and phase information. In a store data step 440, the 
channel data is stored in a channel data buffer. This channel data buffer is located in 

10 memory such as RAM, magnetic media, optical media, or the like. 

[0066] An echo area calculation step 450 includes manipulation of the stored channel 
data using multidimensional de-convolution algorithms. These algorithms are 
mathematical techniques that transform the channel data into multidimensional 
echolocation data. Echo area calculation step 450 can generate the multidimensional 

15 echolocation data without using the transmit lines, receive lines, or scan lines that 
characterize the prior art. 

[0067] Method 400 continues with a store echolocation data step 460 wherein the 
resulting echolocation data are stored using an echolocation data array that utilizes a 
pre-selected coordinate system. The echolocation data is typically located in memory 
20 such as RAM, magnetic media, optical media, or the like. 

[0068] In a step 465, method 400 tests whether the data collection process is complete 
(e. g. the data required to generate the desired image has been collected). If the data 
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collection process is incomplete the method returns to broad-beam design step 410 
wherein another broad-beam is designed. If, at step 465, the data collection process is 
complete an image may be generated in an optional generate image step 470 and 
displayed, on a display device such as a computer monitor, in an optional display step 
5 480. 

[0069] In an alternative embodiment, broad beam design step 410 includes calculation 
of characteristics for several broad-beams. In this embodiment a return to step 410, 
between steps 465 and 420 is optional. The method may proceed directly from step 465 
to transmit step 420 because the desired characteristics for a next broad-beam are pre- 

10 calculated in a prior instance of step 410. 

[0070] FIG. 5 shows a broad-beam system according to an embodiment of the 
invention and generally designated 500. A waveform generator 510, such as a 
programmable pulse sequence generator or the like, is used to generate electronic 
signals, such as electronic pulses 210, that are later used to form a broad-beam 

15 ultrasound beam. The electronic signals are individually delayed, through a delay 

device 515, in several signal channels with a set of delays that are configured to generate 
an ultrasound beam with characteristics designed in step 410 of FIG. 4. The output of 
delay device 515 is coupled to a power amplifier 520, such as a power transistor, 
operational amplifier, high speed FET, or the like, where it is amplified and passed 

20 through a transmit/receive switch 525. Transmit/receive switch 525 optionally includes 
a multiplexer 527 configured to couple input channels including signals received from 
delay device 515 to output channels for transmission to a transducer array 530, which 
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may be analogous to prior art element array 105. Transducer array 530 includes 
ultrasound transducer elements, such as ultrasound transducer elements 110A-110H, 
that generate a broad-beam by converting electrical signals received from 
transmit/receive switch 525 to ultrasound pulses. 
5 [0071] Transducer array 530 is configured to transmit the broad-beam into a material 
under investigation 535. The transmission of the broad-beam occurs in step 420 of FIG. 
4. Echoes are generated in material under investigation 535 through interactions 
between the broad-beam and ultrasound reflective objects, such as tissue and bone. 
Transducer array 530 receives the generated echoes and produces corresponding 
10 electrical signals in step 430 of FIG. 4. These electrical signals, which are typically analog 
electrical signals, are coupled through transmit/receive switch 525 to a variable gain 
amplifier 540, such as a voltage regulated operational amplifier, digitally controlled 
amplifier, amplifying transistor circuit, or the like. 

[0072] After amplification, signals are passed through an optional analog filter 545 to 
15 an A/D converter 550, where the amplified signals are digitized. Analog filter 545 may 
be any analog filter known in the art such as a band-pass filter, a notch filter, or the like. 
A/D converter 550 is typically a commercially available analog to digital converter, or 
the like. 

[0073] The resulting digital data are stored, in step 440 (FIG. 4), in a channel data 
20 storage buffer 555 where they are operated on by signal processor 560. Channel data 
storage buffer 555 may be located in any storage system known in the art. For example, 
channel data storage buffer 555 is optionally located in electronic memory, such as RAM, 
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or magnetic or optical memory such as disc drives, compact disks, or the like. The 
operations performed by signal processor 560 include echo area calculations, of step 450 
(FIG. 4), that transform time domain data stored in channel data storage buffer 555 to 
echolocation data, such as raw data or detected data, that is stored, in step 460 (FIG. 4), 
5 in an echolocation data storage 565. From echolocation data storage 565, data is 
optionally transferred to an additional data storage 570, or accessed by an image 
converter 575. Echolocation data storage 565 and additional data storage 570 may be 
any suitable store devices such as electronic memory, magnetic or optical media, or the 
like. Image converter 575 is analogous to "image scan converters" of the prior art, but 

10 may additionally operate on data generated using a single ultrasound beam rather than 
data generated using a "scan" including several ultrasound beams. In step 470 (FIG. 4), 
image converter 575 may use data stored in echolocation data storage 565, additional 
data storage 570, or both to generate detected data or image data. 
[0074] The image generation process may be analogous to prior art techniques of 

15 image generation using echolocation data generated through beam forming methods. 
For example, a specific position in echolocation data storage 565 is optionally mapped to 
a specific location on a display screen. Intensity and/or color of a position within the 
image may indicate the intensity or other characteristic of echoes detected from within 
material under investigation 535. This image is optionally shown, in step 480 (FIG. 4), 

20 on a display 580 such as an LCD screen, CRT screen, computer monitor, electronic 
display, or the like. 
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[0075] Data used by image converter 575 may result from a series of ultrasound 
beams or alternatively from a single ultrasound beam. Data in additional data storage 
570 is coupled to other components of broad-beam system 500 such as image converter 
575, communications electronics 585 and user interface electronics 590. Components of 
5 broad-beam system 500 are controlled and coordinated by control electronics 595 
through connections not shown in FIG. 5. Control electronics 595 include 
microprocessors, DSPs, and optional computer code 596 configured to control elements 
of broad-beam system 500 and execute methods of the invention such as broad-beam 
process 400. 

10 [0076] FIG. 6 is a flow chart illustrating broad-beam design step 410 according to an 
embodiment of the invention. In this embodiment, calculations are performed using 
computer code 596 and may include, for example, mathematical models of ultrasound 
beam generation, propagation and echoing. In some instances lookup tables are used to 
speed the calculation process. For example, if a user has indicated a specific depth of 

15 analysis a desirable intensity is optionally determined from a lookup table. Broad-beam 
design step 410 begins with a coverage determination step 610 in which the area (or 
volume) within material under investigation 535 to be investigated and the time period 
over which the investigation is to occur is determined. Coverage determination step 610 
may be responsive to options selected by a user and the requirements of the current 

20 imaging (analysis) mode. For example, in a Doppler imaging mode the user may choose 
continuous monitoring and a broad-beam characterized by a continuous series of 
ultrasound pulses. In another example, a user may choose to spotlight a region within 
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material under investigation 535 using a restricted field of view. The choice of a specific 
field of view is optionally used when calculating a width of a generated broad-beam. 
For example, widths of broad-beams may be selected such that an integral number of 
broad-beams fit, with 10% overlap, into a chosen field of view. 
5 [0077] Also, coverage determination step 610 may determine a number of broad- 
beams required to image an area (or volume) within material under investigation 535. 
For example, in one embodiment coverage determination step 610 includes a calculation 
configured to simulate coverage in the far field that determines that an area is best 
imaged using three broad-beams displaced from each other using block-switching 
10 techniques. In other embodiments the calculation determines that an area is best imaged 
using one, two or more broad-beams. When the user has selected a mode of operation 
that includes several different broad-beams, repeated imaging or continuous 
monitoring, coverage determination step 610 is optionally performed once for each 
broad-beam. 

15 [0078] Coverage determination step 610 is followed by a characteristic determination 
step 620 in which further characteristics of broad-beam(s) determined in coverage 
determination step 610 are specified. These characteristics include, but are not limited 
to, ultrasound frequencies, direction, dispersion, pulse shape, phase relationships, 
aperture, intensity, duration, repetition rate and/or other properties of an ultrasound 

20 beam. The characteristics are typically dependent on the imaging mode of analysis 

being performed, the required resolution, and options selected by a user. For example, a 
continuous monitoring mode may require a broad-beam generated at a specific pulse 
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rate, high resolution may require use of multiple ultrasound frequencies, and a user may 
choose to investigate a narrow region best probed by a broad-beam with low dispersion. 
In addition to the characteristics discussed above, characteristic determination step 620 
may include selection of a coordinate system with which to represent the area covered 
5 by the broad-beam and an origin of this coordinate system. Such a coordinate system 
may be used to store echolocation data. Selection of a coordinate system is optionally 
responsive to the shape of a broad-beam. Examples of possible coordinate systems are 
illustrated in FIG. 7. 

[0079] Coverage determination step 610 and characteristic determination step 620 are 
10 optionally responsive to resolution and dynamic range requirements. For example, in 
one embodiment these steps are responsive to user input that specifies an image zoomed 
in on a specific area. In another embodiment these steps are responsive to user input 
that specifies a higher image resolution for part or all of an image. In another 
embodiment coverage determination step 610 includes a determination that a single 
15 ultrasound beam should be generated but that, for instance to enhance resolution, the 
echoes generated by the single ultrasound beam should be detected by several different 
sets of receive transducers in multiple transmit/receive cycles. 

[0080] Coverage determination step 610 and characteristic determination step 620 are 
optionally responsive to feedback generated in other steps of the invention. For 
20 example, in one embodiment, echolocation data indicates that a region of the covered 
area is poorly imaged and that the poor imaging is caused by a highly reflective 
boundary disposed between the poorly imaged region and the closest of transducer 
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elements 110. In response to this feedback, coverage determination step 610 and 
characteristic determination step 620 include defining a steered broad-beam that probes 
the region from alternative ultrasound transducers that are not inline with the reflective 
boundary and the region to be probed. 
5 [0081] Broad-beam selection step 630 includes selection of a broad-beam for 
transmission. The broad-beam is selected from those defined in characteristic 
determination step 620. If several broad-beams have been characterized in characteristic 
determination step 620 then broad-beam selection step 630 is optionally performed more 
than once before the next occurrence of characteristic determination step 620. In such a 

10 case broad-beam selection step 630 is repeated after step 465 of FIG. 4. 

[0082] Broad-beam design step 410 concludes with a calculate excitation step 640. 
Calculate excitation step 640 includes determining the proper physical parameters 
required to generate the broad-beam selected in broad-beam selection step 630. These 
physical parameters include, for example, which transducer elements 110 to excite, 

15 electronic pulse voltages, pulse delay times, multiplexer 527 settings, and/or the like. 
For example, in one embodiment a selected ultrasound beam, having a particular 
desired shape and direction, requires use of a specific set of transducer elements 110, 
excited by a particular electronic waveform characterized by amplitudes, frequencies 
and phases, each of the required set of transducer elements 110 being excited with an 

20 appropriate delay. The proper physical parameters are determined, for example, using 
a mathematical model to calculate a voltage, waveform, and delay used for exciting a 
particular member of transducer elements 110. In one embodiment the voltage is 
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responsive to a distance into the material under investigation 535 the broad-beam is 
expected to penetrate. 

[0083] FIGS. 7A-7C show embodiments (710A - 710C) of a broad-beam 710 
determined in coverage determination step 610 and characteristic determination step 
5 620. FIG. 7 A shows broad-beam 71 OA generated using a linear embodiment of 

transducer array 530. The area of an insonified region, generally designated 715A, is 
optionally represented by a radial (0,R) coordinate system with an origin 720 located at 
the surface of transducer elements 110. Points within insonified region 715 are identified 
by their distance (R) from an origin 720 and their angular coordinate (9) relative to 

10 transducer array 530 or an axis, such as an axis 730 or an axis 735. In alternative 

embodiments the focal point of broad-beam 710B is located behind transducer array 530, 
rather than in front of transducer elements 110 as shown in FIG. 2. 
[0084] FIG. 7B shows broad-beam 710B generated using a curvilinear embodiment of 
transducer array 530. An insonified region, generally designated 715B, is optionally 

15 represented by a radial coordinate system with an origin 755 behind transducer array 
530. This origin location provides insonification of more area proximal to transducer 
elements 110 than an origin location closer to transducer array 530 as shown in FIG. 7A. 
The location of origin 755 behind transducer array 755 is optionally independent of the 
shape of transducer array 755. Embodiments of the invention also include, but are not 

20 limited to, positioning origin 755 and/or a focal point behind a linear embodiment of 
transducer array 530. 
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[0085] FIG. 7C shows broad-beam 71 0C that results in an insonified region, generally 
designated 715C. Insonified region 715C is more rectangular in shape than those 
generated by broad-beam 71 OA and broad-beam 71 OB, shown in FIGS. 7 A and 7B, 
respectively. The region insonified by broad-beam 710C may be preferably represented 
5 by a Cartesian (x,y) coordinate system 780 because of the region's rectangular shape. 
[0086] In contrast with the prior art, where the maximum intensity is found at the 
center of an ultrasound beam, the maximum intensity of a broad-beam, such as broad 
beam 710B or 710C, may be at points other than along the beam's center. FIG. 7D shows 
a plot 790 of ultrasound intensity through a cross-section of broad-beam 71 0C as 

10 measured at a distance from transducer array 530, approximately equal to 1/2 the width 
of the beam's aperture. This cross-section is indicated by a dashed line 785 in FIG. 7C. 
In some circumstance, the intensity profile of a broad-beam represents a more desirable 
energy distribution than those found in the prior art. For example, the energy 
distribution illustrated by plot 790 is more evenly distributed over insonified region 

15 71 5C than the energy distribution within a prior art ultrasound beam in the region of a 
focal point. 

[0087] FIG. 8 shows details of an embodiment of transmit step 420 of FIG. 4. In this 
embodiment, step 420 includes a waveform generation step 810 in which waveform 
generator 510 is used to generate an electrical waveform with characteristics calculated 
20 in broad-beam design step 410. The generated waveform optionally includes a plurality 
of pulses of varying frequency or phase. In a signal delay step 820 the generated 
waveform is reproduced in several signal channels and delayed, using delay device 515, 
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by times determined in broad-beam design step 410. Waveforms in each signal channel 
are amplified in an amplification step 830 using power amplifier 520. The amplified 
waveforms are coupled through multiplexer 527 in a multiplex step 840. Multiplexer 
527 is set to direct the waveform in each signal channel to one or more member of 
5 transducer elements 110 in transducer array 530. In sound generation step 850, the 
directed waveforms cause transducer array 530 to generate broad-beam 710, which is 
directed into material under investigation 535. Sound generation step 850 completes 
transmit step 420. 

[0088] FIG. 9 shows details of an embodiment of receive step 430 of FIG. 4 in which 
10 echoes are detected and converted to digital data. In a set switch step 910 

transmit/receive switch 525 is set such that signals produced at transducer elements 110 
are coupled through multiplexer 527 to variable gain amplifier 540. In an echo detection 
step 920, echoes from within material under investigation 535 are detected by members 
of transducer elements 110 in transducer array 530. The members of transducer 
15 elements 110 used for detection of echoes are optionally different than the members of 
transducer elements 110 used to transmit broad-beam 710. In various embodiments 
these two sets of transducer elements 110 are configured a number of ways. For 
example the sets may be identical, interleaved, overlapped partially along transducer 
array 530 or not overlapped along transducer array 530. The electronic signals resulting 
20 from the detected echoes are coupled to variable gain amplifier 540 because 
transmit/receive switch 525 was set in set switch step 910. 
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[0089] The electronic signals coupled to variable gain amplifier 540 are amplified in a 
variable amplification step 930. Variable amplification step 930 optionally includes 
feedback based on data obtained using a prior broad-beam 710. The feedback provides 
adaptive processing and can be used to adjust signal within each channel such that the 
5 dynamic range of subsequent data manipulation steps are maximized. For example, in 
one embodiment, if previous execution of variable amplification step 930 resulted in the 
saturation of a specific channel, then amplification in that channel is optionally reduced 
in a following execution of variable amplification step 930. The reduction, or adaptive 
front end gain, is compensated for in later data manipulation that occurs after 
10 digitization of the amplified signal. In another embodiment, transducer elements 110 
near the center of transducer array 530 are found to systematically respond to echoes 
more strongly than transducers elements 110 near an edge of transducer array 530. 
Variable amplification step 930 optionally includes compensation for this systematic 
difference. 

15 [0090] In an optional analog filtering step 940 the electronic signals, amplified in 
variable amplification step 930, are processed using analog filter 545. This processing 
includes, for example, I/Q mixing, removal of unwanted frequencies and shifting of 
signals into frequency ranges more suitable for further data manipulation. 
[0091] In a data conversion step 950 the electronic signals, optionally filtered in analog 

20 filtering step 940, are digitized using A/D converter 550. The generation of digital data 
completes receive step 430 (FIG. 4). In various embodiments data conversion step 950 
occurs at alternative times within broad beam process 400. After the completion of 
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receive step 430 the resulting digital data is stored, in store data step 440 (FIG. 4), in 
channel data storage buffer 555. 

[0092] FIG. 10 shows an embodiment of a channel data array 1000 configured to hold 
the digital data stored in store data step 440. Channel data array 1000 is stored in 
5 channel data storage buffer 555. A first axis 1010, of Channel data array 1000, is indexed 
by echo receiving members of transducer array 530. A second axis 1020 of channel data 
array 1000 is divided into time channels. Values stored at each location in the array 
indicate the intensity and phase of echo signals detected by a specific member of 
transducer array 530 at a specific time. 

10 [0093] Channel data storage buffer 555 optionally includes several channel data array 
1000. Additionally, the information stored in channel data array 1000 may be used to 
average or sum received signals. In various embodiments channel data array 1000 is 
configured to store multidimensional data. For example, in one embodiment transducer 
array 530 is a two dimensional array of transducer elements 110. In this embodiment 

15 channel data array 1000 includes two axes representing the two dimensions of 
transducer array 530 and one axis representing time channels. 
[0094] Echo area calculation step 450 uses data stored in store data step 440 to 
generate echolocation data indicating the positions and strengths of echo sources within 
material under investigation 535. This generation of echolocation data includes 

20 transformation of multidimensional time-channel data, within channel data array 1000, 
to multidimensional positional (echolocation) data. For example, in one embodiment 
two-dimensional time-channel data is transformed into echolocation data represented 
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by two-dimensional spatial coordinates. The data transform of echo area calculation 
step 450 is performed using a variety of alternative transform algorithms, examples of 
which are disclosed herein. These transforms are optionally used to generate two- 
dimensional echolocation data using signals received as the result of a single broad- 
5 beam 710. In an alternative embodiment echo area calculation step 450 is replace by an 
analogous echo volume calculation step including an additional spatial dimension. 
Echo Volume calculation includes the generation of three-dimensional echolocation data 
using signals received as the result of a single broad-beam, the broad-beam covering a 
three dimensional volume. 

10 [0095] FIGs. 11A and 11B show two embodiments of an echolocation data array 1100 
stored in echolocation data storage 565 and configured to store positional data resulting 
from echo area calculation step 450. These two embodiments employ different 
coordinate systems. As discussed in further detail below, the more efficient coordinate 
system may be dependent on, among other factors, the shape of an individual 

15 ultrasound beam 710. In most instances, a more efficient coordinate system will overlay 
closely with the area being insonified. For example, as shown in FIGs. 7A-7C, the area 
insonified by broad-beam 71 OA, broad-beam 710B and broad-beam 710C are each 
preferably represented by different coordinate systems with different origins. Use of a 
more efficient coordinate system may increase sampling efficiency and spatial 

20 resolution. Selection of a preferred coordinate system and echolocation data array 1100 
may be responsive to the shape of an ultrasound beam, such as broad-beam 710, and 
optionally occurs in steps 410, 440 or 450. 
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[0096] FIG. 11 A shows an embodiment of echolocation data array 1100 using a 
Cartesian coordinate system including a first axis 1110 indicating an X coordinate 
(position) and a second axis 1120 indicating a Y coordinate (position). FIG. 11B shows 
an alternative embodiment of echolocation data array 1100 using a radial coordinate 
5 system including first axis 1110 indicating an angle (9) coordinate and second axis 1120 
indicating a radius coordinate. Alternative embodiments of echolocation data array 
1100 are represented by alternative coordinate systems. Additional data, not shown, is 
optionally used to relate first axis 1110 and second axis 1120 to transducer array 530. For 
example, echolocation data array 1100 is optionally characterized by vectors relating the 
10 origin of each coordinate system to a specific member of ultrasound transducer elements 
110. 

[0097] FIGs. 12 A and 12B illustrate how use of one coordinate system may be more 
efficient than use of another coordinate system. FIGs. 12A and 12B show the 
embodiments of echolocation data array 1100 shown in FIGs. 11 A and 11B, respectively, 

15 overlaid on an ultrasound beam 1210. Ultrasound beam 1210 is an embodiment of 
broad-beam 710. FIG. 12A shows a Cartesian coordinates system including, for the 
purposes of illustration, eleven "X" divisions separating data bins 1220. Data bins 1220 
are just adequate to cover the far field, generally designated 1230. Because the spacing 
of data bins 1220 in the X dimension is the same in the near field, generally designated 

20 1240, a number of data bins 1220 in near field 1240 are mapped to area that is not probed 
by ultrasound beam 1210. These data bins 1220, not mapped to probed area, represent 
inefficient sampling of the material under investigation 535. 
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[0098] In contrast, FIG. 12B shows use of a radial coordinates system to represent the 
area insonified by ultrasound beam 1210. In the radial coordinate system the size of 
data bins 1250 vary as a function of the "R" coordinate. Data points in this embodiment 
of echolocation data array 1100 are, therefore, more efficiently mapped to the area 
5 probed by ultrasound beam 1210, than the embodiment of echolocation data array 1100 
shown in FIG. 12A. The variation of data bin 1250 size increases efficiency because, as 
shown in FIG. 12B, a greater fraction of data bins 1250 within data array 1100 fall within 
the area covered by ultrasound beam 1210. 

[0099] The granularity of data bins is dynamic. In some embodiments echolocation 
10 data array 1100 represents a Nyquist sampled space wherein the density of bins 1250 is 
varied such that the number of samples just satisfies Nyquist criteria for un-aliased 
sampling throughout a region of interest. In some embodiments the density of bins 1250 
is varied such that the resolution of resulting echolocation data is greater in a specific 
region. For example, in one embodiment a user specifies a particular region where more 
15 image detail is desired. In response, broad-beams systems and methods use an 
echolocation data array 1100 with greater density of bins 1250 in this region. 
[00100] Some embodiments of the present invention include extrapolation and 
interpolation between data bins 1250. For example, in one embodiment interpolation is 
used in the far field, where each of data bins 1250 represent a greater area, to increase 
20 the density of echolocation data. Optionally, less interpolation is used in the near field 
were the density of data bins 1250 is greater. 
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[00101] The resolution (sampling frequency) of channel data generated in receive step 
430 fundamentally limits the resolution of resulting echolocation data as a result of the 
Nyquist theorem. However, the resolution of data generated in receive step 430 is 
optionally improved through signal averaging or up-sampling techniques. Up- 
5 sampling techniques include the use of additional data and optionally include feedback 
such that additional data is collected in regions where improved resolution is most 
needed. 

[00102] FIGs. 13 through 15 are used to show embodiments of echo area calculation 
step 450 (FIG. 4). FIG. 13 shows propagation of ultrasound between transducer elements 

10 110A-110S, and ultrasound reflecting objects within material under investigation 535. 
FIG. 14 shows channel data produced from detected echoes. And, FIG. 15 shows 
echolocation data generated using the channel data shown in FIG. 14. 
[00103] In several embodiments of echo area calculation step 450 including data 
transform methods it is assumed that the primary contributor to detected echoes from 

15 each location within the material under investigation 535 is the member of transducer 
elements 110 closest to that location. This element is referred to as the main contributing 
element (MCE). Typically, the member of transducer elements 110 that is closest to a 
location is the MCE for that particular location, and any ultrasound reflective object at 
that location. However, the identity of the MCE may also be dependant on the direction 

20 of broad-beam 710 and the shape of transducer array 530. In such a case, the MCE may 
not be the transducer element 110 closest to the particular location. The data transform 
methods, of echo area calculation step 450 (FIG. 4),optionally include broad-beam 710 
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direction, transducer array 530 geometry, feedback, as well as other factors for 
determining an MCE that is not the closest member of transducer elements 110 to an 
ultrasound reflective object. 

[00104] FIG. 13A shows ultrasound 1305 transmitted from a single transducer element 
5 HOG. Ultrasound 1305 travels through material under investigation 535 (not shown) 
until it strikes an ultrasound reflecting object 1310A. Transducer element HOG is the 
closest of transducer elements 110A-110S to ultrasound reflecting object 1310A, and is 
therefore considered to be the MCE for reflecting object 1310A. At ultrasound reflecting 
object 1310A, ultrasound 1305 generates ultrasound echoes 1315 of which ultrasound 
10 echoes 1315A-1315F are shown. Ultrasound echoes 1315 propagate back to transducer 
elements 110A-110S where they are detected. 

[00105] Although FIG. 13A shows ultrasound 1305 transmitted from one transducer 
element HOG (the MCE), in most embodiments ultrasound is transmitted from a 
plurality of transducer elements 110A-110S during the formation of broad-beam 710. 
15 FIG. 13B shows ultrasound 1330 generated by a single transducer element 110Q, which 
is the MCE for an ultrasound reflecting object 1310B. Echoes 1340, of which ultrasound 
echoes 1340A-1340F are shown, generated at reflecting object 1310B travel back to and 
are detected by transducer elements 110A-110S. 

[00106] FIG. 14 shows an embodiment of channel data array 1000 including data 
20 generated by ultrasound 1305 and ultrasound 1330 shown in FIG. 13. Each of columns 
1410A-1410S in channel data array 1000 represents signal(s) detected at one of 
transducer elements 110A-110S. Each of rows 1420A-1420U in channel data array 1000 



includes the signal detected during a specific time period. In FIG. 14 data elements 1430, 
that included data generated by detection of echoes 1315 and 1340, are those data 
elements 1430 that intersect a data location line 1440A or a data location line 1440B, 
respectively. Thus, ultrasound echoes generated from a reflective object, such as 
5 ultrasound reflective object 1310, within material under investigation 535 results in data 
that lies along a line, such as data location lines 1440A or 1440B. Data location lines 
1440A and 1440B can be calculated from first principles of physics and geometry using a 
known geometry of transducer array 530 and the speed of sound within material under 
investigation 535. Data location lines 1440A and 1440B do not intersect the MCE, 

10 transducer element HOG, nor typically any other transducer element 110. In practice, 
material under investigation 535 includes numerous ultrasound reflective objects 1310, 
and channel data array 1000 includes data generated by each. 
[00107] In embodiments of echo area calculation step 450, echolocation data is 
calculated by summing data along a line such as data location line 1440A, data location 

15 line 1440B, or the like. For example, summation of data along data location line 1440B 
generates a result indicative of the magnitude of echoes 1315 generated at the position 
occupied by ultrasound reflecting object 1310B and represented by a data bin, such as 
data bin 1220 or data bin 1240. The sum is stored in the representative data bin. A 
similar summation is optionally performed for each data bin in echolocation data array 

20 1100. Through multiple summations echolocation data array 1100 is populated with 
echolocation data representing ultrasound reflective objects within material under 
investigation 535. 
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[00108] FIG. 15 shows an embodiment of echolocation data array 1100 including 
echolocation data bins 1520. Each of echolocation data bins 1520 is associated with a 
unique line, such as data location line 1440A, in channel data array 1000 as shown in 
FIG. 14. Data along the unique line is summed to calculate the magnitude of echo 
5 generation that occurred at the physical locations represented by each of data bins 1520. 
This summation is optionally performed for all of data bins 1520 and thus can be used to 
calculate echolocation data over the entire echolocation data array 1100. 
[00109] FIG. 16 shows a data transform method included in an embodiment of echo 
area calculation step 450. This embodiment includes a select element step 1610 in which 

10 one of echolocation data bins 1520, within echolocation data array 1100, is selected. 

Typically, selection of each of echolocation data bins 1520 is accomplished by traversing 
echolocation data array 1100 in a systematic fashion. Select element step 1610 is 
followed by a determine line step 1620 in which the unique line in channel data array 
1100 associated with the selected echolocation data bins 1520 is determined. 

15 Determination is accomplished by calculating the line from geometric principles, using a 
look-up table with previously calculated lines, or the like. Determination may occur 
before or during echo area calculation step 450. In various embodiments determination 
occurs prior to or during broad-beam design step 410. In alternative embodiments, 
determination occurs during steps 420, 430, and/or 440 (FIG. 4). Determine line step 

20 1620 is followed by a sum data step 1630 that includes summation of data from data 
elements 1430 that intersect the line determined in determine line step 1620. In one 
embodiment sum data step 1630 includes a simple addition of data. In alternative 
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embodiments sum data step 1630 includes use of weighting functions, matrix 
manipulation, extrapolation, interpolation, or like mathematical techniques. In one 
embodiment sum data step 1630 is facilitated by firmware within control electronics 595. 
In a store result step 1640 the result of the summation of step 1630 is stored in the data 
5 element selected in select element step 1610. 

[00110] Steps 1610 through 1640 are optionally repeated for all echolocation data bins 
1520 in echolocation data array 1100. FIG. 15 shows two sets (1550 and 1560) of 
echolocation data bins 1520 including non-zero values resulting from summation along 
data location lines 1440A and 1440B using the method shown in FIG. 16. Each set (1550 
10 and 1560) of echolocation data bins 1520 typically include echolocation data bins 1520 
with differing non-zero values. In several embodiments one or more of steps 1610 
through 1640 are performed as parallel processes. 

[00111] Alternative embodiments of echo area calculation step 450 include alternative 
methods of data transformation. These methods use, for example, calculations 

15 performed in the frequency domain, use of phase relationships between received 
signals, use of apodization functions to weigh contributions of each of transducer 
elements 110, feedback mechanisms, correlation analysis and consideration of 
transmitting transducer elements 110 other than the MCE. These other transducer 
elements 110 are used to improve both the quality and speed of the transformation from 

20 channel data to echolocation data. 

[00112] In one embodiment, echo area calculation step 450 includes use of an 
apodization function to weigh contributions of each transducer element 110. Weighting 
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may be desirable because those transducer elements 110 closer to an MCE receive 
stronger echoes from a particular reflective object 1310 than do transducer elements 110 
further from the MCE. Signals detected at an MCE and the transducer elements 110 
nearby are therefore given greater weight than transducer elements 110 further from the 
5 MCE. 

[00113] FIG. 17 shows three alternative apodization functions according to 
embodiments of the invention. Graph 1710 illustrates these three alternative 
apodization functions, designated 1720, 1730 and 1740. For example, if transducer 
element HOG is the MCE for one of data elements 1430 selected in select element step 
10 1610 of FIG. 16, then apodization function 1720 is used in sum data step 1630 such that 
the resulting sum includes a greater contribution from transducer elements 110 near 
transducer element HOG. Likewise, for summations wherein transducer elements 110K 
and 11 OS are the MCE, apodization functions represented by lines 1730 and 1740 are 
optionally used. 

15 [00114] In alternative embodiments, echo area calculation step 450 is performed at least 
in part in the frequency domain. Data is converted using invertible transforms, for 
example sine transform, Fourier transform, wavelet transform, or the like. 
[00115] In some embodiments of echo area calculation step 450 phase relationships 
between received signals are used to distinguish between those signals resulting from 

20 ultrasound transmitted by the MCE and those signals resulting from secondary 

contributing elements (SCEs). SCEs are transducer elements 110, other than the MCE, 
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that contribute to signal arising from a given ultrasound reflective object, such as 
ultrasound reflective object 1310. 

[00116] FIG. 18 shows ultrasound 1810 and 1305 transmitted from transducer elements 
110F and HOG and striking ultrasound reflective object 1310A. Transducer element 
5 HOG is considered the MCE for ultrasound reflective object 1310 A because it is the 
closest member of transducer elements 110. In alternative embodiments a closely 
grouped set of transducers are treated jointly as an MCE. Other transducer elements 
110, such as transducer element 110F, also generate ultrasound that can reach reflective 
object 1310A. In this example, transducer element 110F is a SCE. Ultrasound must 
10 travel further from these (SCE) transducer elements 110 than from the MCE transducer 
elements 110, before reaching ultrasound reflective object 1310A. As with the 
ultrasound generated by the MCE, ultrasound from the SCEs generates echoes when 
striking ultrasound reflective object 1310A. Some of these echoes are detected at 
transducer array 530. 

15 [00117] FIG. 19 shows locations of signals generated by SCE transducer element 110F in 
channel data array 1000. These signals lay along a data location line 1910 similar to data 
location line 1440A, but at a slightly later time. The time difference between data 
location lines 1440A and 1910 is due to the difference in time required for ultrasound to 
travel to ultrasound reflective object 1310A from transducer element 110F and from 

20 transducer element HOG. It is desirable to distinguish data resulting from SCEs from 
data resulting from an MCE. Although signal from the MCE is typically stronger than 
signal resulting from SCEs (due to the longer distance ultrasound must travel), the 
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signal from the SCEs is additionally differentiated by a phase difference that results 
from the difference in distance traveled. Considering signals only with specific phases 
allows signals resulting from SCEs to be separated by filtering. For example, in one 
embodiment SCE signal is filtered out by more than 10 dB and in some embodiments by 
5 more than 38 dB. 

[00118] In various embodiments, data resulting from SCEs are used to improve results 
obtained in echo area calculation step 450. For example, in some embodiments, data 
resulting from an SCE is added to data resulting from an MCE. Thus, data laying along 
data location line 1910, as shown in FIG. 19, is added to data lying along data location 

10 line 1440A. The data lying along data location line 1910 includes data resulting from 
ultrasound generated at (SCE) transducer element 110F and echoed from reflecting 
object 1310A. After a phase adjustment and weighting this data may constructively add 
to data lying along data location line 1440A, and thus improve the signal to noise ratio of 
echolocation data indicating the presence of reflecting object 1310A. Typically, SCEs 

15 closest to an MCE are given more weight than SCEs further away. For example, one 
embodiment uses a Chi Squared weighting distribution, centered on the MCE to 
determine weighting of neighboring SCEs. In another embodiment the weighting 
distribution is responsive to feedback algorithms that reduce the weight of SCEs whose 
signal in channel data array 1000 overlap with a strong MCE signal. 

20 [00119] In other embodiments signal resulting from an SCE is subtracted from signal 
resulting from an MCE. For example, if a large MCE signal is detected along data 
location line 1440A as shown in FIG. 19, then a correspondingly large SCE signal will be 
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expected along data location line 1910. Since the corresponding SCE signal is 
predictable and approximate values can be calculated as a function of the MCE signal, 
the calculated values can be subtracted from channel data values stored in data elements 
1430 before these data values are used to calculate values for other echolocation data 
5 bins 1520. Consideration of data resulting from SCEs to improve echo area calculations 
optionally occur as part of sum data step 1630 (FIG. 16). 

[00120] Several embodiments of echo area calculation step 450 use feedback. For 
example, in one embodiment calculated echolocation data is processed in a "reverse 7 ' 
data transform using techniques that produce a simulated echo signal (simulated 

10 channel data) based on the calculated echolocation data. This reverse transform 

produces a simulation of the channel data that would be expected if the calculation of 
echolocation data was optimal. The reverse transform is optionally preformed using 
ray-tracing methods known in the art. The simulated channel data is compared with the 
actual echo data stored in channel data array 1000. Similarity between these two data 

15 sets is indicative of the quality of the calculation used to produce the echolocation data. 
In an optional iterative process, the calculation of echolocation data is repeated using 
varying parameters responsive to this similarity. These parameters may include 
different weighting factors, apodization functions or SCEs, manipulated to optimize the 
similarity between the data in channel data array 1000 and simulated echo signals. 

20 [00121] In other embodiments feedback includes use of echolocation data to control 
broad beam design step 410. For example, in one embodiment the direction of an 
ultrasound beam designed in step 410 is responsive to the location of reflective 
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boundaries in material under investigation 535. In other examples, the focus, width, 
frequency, intensity, or number of beams designed in step 410 are responsive to 
calculated echolocation data. 

[00122] Several embodiments of echo area calculation step 450 include data transforms 
5 employing correlation analysis. Correlation methods are known in the data analysis art 
and are useful for enhancing similarities and making comparisons between data. 
Correlation is particularly useful for comparing data that systematically differs, for 
example by a change in phase. A cross-correlation analysis of two data sets, differing by 
a constant degree along one coordinate, identifies the constant difference and the 

10 similarity of the data after accounting for that difference. An auto-correlation analysis of 
a data set exemplifies periodic or repetitive signals within the data. 
[00123] FIG. 20 shows an embodiment of echo area calculation step 450 that includes a 
cross-correlation method used to identify components of SCE data that correlate well 
with MCE data. In a calculate cross-correlation step 2010 data laying along a line, such 

15 as line data location 1440 A (FIG. 14), associated with an MCE is cross-correlated with 

data laying along a line, such as data location line 1910 (FIG. 19), associated with an SCE. 
Each of these sets of data is optionally pre-processed using a function such as 
apodization function 1720. The cross-correlation generates a correlation data set that can 
be expressed as a function of phase difference verses similarity between the two data 

20 sets. In a calculate phase difference step 2020 the expected phase difference between the 
MCE data and the SCE data is calculated based on a known geometrical relationship 
between the MCE and the SCE. In a look-up step 2030 this calculated phase difference is 



used to look-up a similarity value in the correlation data set generated by the cross- 
correlation, at that specific phase difference, in the correlation data set. The similarity 
value, corresponding to the phase difference calculated in phase difference step 2020, is 
indicative of how useful the SCE data can be in improving the signal to noise ratio of the 
5 MCE data because more similar SCE data is more likely to constructively add to the 
MCE data. In a decision step 2040 the similarity value is compared with a 
predetermined threshold. If the similarity value is greater than the threshold then the 
SCE data is added to the MCE data in an add data step 2050. If, in step 2040, the 
similarity value is found to be less than the predetermined threshold, computer code 596 

10 determines, in a decision step 2060, if further analysis of the particular SCE data set is 

warranted. Further analysis may be warranted if, for example, nearby SCEs are yet to be 
examined or if a user has requested additional improvement in the signal to noise ratio. 
If not, the analysis of this particular SCE data set is concluded. If step 2060 determines 
that further analysis is warranted then the SCE data set is processed in an optional filter 

15 step 2070. The processing in step 2070 includes filtering, truncation or similar means 
designed to enhance the components of the SCE data set that correlate well with the 
MCE data set. For example, in one embodiment an alternative function, such as 
apodization function 1740 is applied to the SCE data set. The steps shown in FIG. 20 are 
optionally applied to more than one SCE data set. 

20 [00124] Echolocation data generated using alternative embodiments of echo area 
calculation step 450 are optionally compared, and the comparison may be used to 
determine the consistency of calculations or to provide feedback. For example, in one 
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embodiment two repetitions of echo area calculation step 450 include consideration of 
different SCEs. The accuracy of these calculations is checked by comparing the results of 
each repetition. The closer the results the more likely the use of SCEs is producing an 
accurate result. In another example, echolocation data calculated using two different 
5 embodiments of echo area calculation step 450 are found to be significantly different. 
These differences are used as feedback affecting other steps in the broad-beam 
technology. For example, ^reproducibility of echolocation data in a specific region is 
optionally used to provide feedback to broad-beam design step 410 such that a 
characteristic (intensity, frequency, direction, etc.) of a broad-beam probing that region 
10 is modified. 

[00125] Data stored in echolocation data array 1100 is optionally used in generate 
image step 470 (FIG. 4) generate images of material under investigation 535 that can be 
displayed to a user. This generation and display is accomplished using image converter 
575 and display 580, respectively. Since two dimensional data can be generated from a 

15 single ultrasound beam using broad-beam techniques a two dimensional image can be 
generated from a single ultrasound beam. In various embodiments this capability 
increases the image frame rate relative to prior art methods because an image is 
produced in a time limited by a single pulse return time, or optionally the return time of 
a few pulses (i.e. <5, <10, <20, <40 or <64), rather than the many (i.e. >100) pulse return 

20 times of the prior art. Benefits of generating an image from a single ultrasound beam 
include possibly reducing jitter in the resulting image because, relative to the prior art, 
there is less time for relative movement between transducer array 530 and material 

49 



under investigation 535 during the period data is collected. Generating an image from a 
single ultrasound beam may also reduce the amount of ultrasound energy directed into 
material under investigation 535 and the amount of electrical power required to generate 
an image. 

5 [00126] From the description of the various embodiments of the process and apparatus 
set forth herein, it will be apparent to one of ordinary skill in the art that variations and 
additions to the embodiments can be made without departing from the principles of the 
present invention. For example, transducer elements 110 can be replaced by alternative 
ultrasound generating elements and transmit/receive switch 515 can be replaced by 
10 separate transmit and receive switches. The number of transducer elements 110 shown 
in the figures is not meant to be limiting. Typical embodiments include larger numbers 
of transducer elements 110. Likewise, the resolution of shown data arrays is selected for 
illustrative purposes only. Typical embodiments include data arrays with larger 
numbers of data bins. 

15 [00127] Broad-beam technology is applicable to systems configured to use both area 
forming and conventional beam forming. Some embodiments include means for 
switching between theses two approaches. For example, area forming may be used to 
survey and area and conventional beam forming techniques may be used to focus 
energy onto a specific area of interest. In some embodiments, including two 

20 dimensional transducer arrays, area forming is used at the same time as conventional 
beam forming techniques. For example, one set of transducer elements may be used for 
area forming while another set of transducer elements may be used for conventional 
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beam forming. In another example, area forming may be used to gather data in one 
spatial dimension while conventional beam forming is used to gather data in an other 
spatial dimension. An ultrasound beam may be configured for area forming in one 
dimension and conventional beam forming in another dimension. In these examples, 
more than one method of echolocation is performed at the same time, each method 
optionally being associated with a specific spatial dimension. 
[00128] Broad-beam technology is applicable to any system limited by the use of 
phased arrays to scan a focused beam over an area or volume. These systems may 
include sonic systems such as sonar, as well as electromagnetic systems such as radar. 
Embodiments of broad-beam technology are used with two dimensional transducer 
arrays. In these embodiments echo volume calculations replace echo area calculations 
and the transform of step 450 includes conversion of a three dimensional (Transducer, 
Transducer, Time) array of echo data to a three dimensional (x, y, z) echolocation data 
array. In one embodiment a single three dimensional ultrasound beam is used to 
perform volume forming and thus produce echolocation data covering a volume in 
space. 
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